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ABSTRACT 

We have analyzed the full set of SOHO/MDI f- and p-mode oscillation frequencies 
from 1996 to date in a search for evidence of solar radius evolution during the rising 
phase of the current activity cycle. Like Antia et al. (2000), we find that a significant 
fraction of the f-mode frequency changes scale with frequency; and that if these are 
interpreted in terms of a radius change, it implies a shrinking sun. Our inferred rate of 
shrinkage is about 1.5 km/y, which is somewhat smaller than found by Antia et al. We 
argue that this rate does not refer to the surface, but rather to a layer extending roughly 
from 4 to 8 Mm beneath the visible surface. The rate of shrinking may be accounted 
for by an increasing radial component of the rms random magnetic field at a rate that 
depends on its radial distribution. If it were uniform, the required field would be ~ 7 
kG. However, if it were inwardly increasing, then a 1 kG field at 8 Mm would suffice. 

To assess contribution to the solar radius change arising above 4Mm, we analyzed the 
p-mode data. The evolution of the p-mode frequencies may be explained by a magnetic 
field growing with activity. Our finding here is very similar to that of Goldreich et al. 
(1991). If the change were isotropic, then a 0.2 kG increase, from activity minimum 
to maximum, is required at the photosphere, which would grow to about 1 kG at 1 
Mm. If only the radial component of the field were to increase, then the requirement 
for the photospheric field increase is reduced to a modest 60-90 G. A relative decrease 
in temperature of the order of 10 -3 in the sub-photospheric layers, or an equivalent 
decrease in the turbulent energy, would have a similar effect to the required inward 
growth of magnetic field change. 

The implications of the near-surface magnetic field changes depend on the anisotropy 
of the random magnetic field. If the field change is predominantly radial, then we infer 
an additional shrinking at a rate between 1.1-1.3 km/y at the photosphere. If on the 
other hand the increase is isotropic, we find a competing expansion at a rate of 2.3 km/y. 
In any case, variations in the sun's radius in the activity cycle are at the level of 10 -5 
or less, hence have a negligible contribution to the irradiance variations. 



Subject headings: Sun: radius — Sun: activity — Sun: oscillations — Sun: interior 
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1. Introduction 

Measuring the sun's radius, and its variability, are 
significant, long-standing problems, especially in the 
context of understanding the cause of solar irradiance 
variations. Recently, it has been pointed out that he- 
lioseismology can provide a useful measure of the solar 
radius. Schou ct al. (1997) and Antia (1998) showed 
that f-modes frequencies are good probes of the ra- 
dius, and they inferred a value of the solar radius 
which is about 300 km smaller than the one adopted 
in solar models at that time. The model values were 
based on a direct measurement of the sun's photo- 
spheric radius. The smaller radius has been con- 
firmed by Brown and Christensen-Dalsgaard (1998) 
from many years of transit measurements using the 
Solar Diameter Monitor. The connection between the 
"true" solar radius and that inferred from f-modes is 
explained in Section 3.1. 

Following the suggestion of Schou et al. (1997), 
Dzicmbowski ct al. (1998, 2000) used f-modc data 
from the MDI/SOHO instrument to determine the 
evolution of the seismic solar radius through the ris- 
ing phase of the present activity cycle. They reported 
statistically significant variations that showed no ap- 
parent correlation with activity measures. On the 
other hand, with GONG f-mode frequencies (rang- 
ing between 1.015 mHz to 1.425 mHz, or equivalently 
£ from 100 to 200), covering the rising phase of ac- 
tivity to the beginning of 1999, Antia et al. (2000) 
found a net decrease of about 5 km in the solar radius. 
They attributed the difference with Dziembowski et 
al. (1998) to the latter's use of higher degree modes 
(up to ^=300). They pointed out the latter ^-modes 
are more likely to be effected by factors other than an 
evolving radius. 

In the present work, we use oscillation data from 
SOHO/MDI covering 1996.3-2000.5. We first use f- 
mode data in an effort to infer a signal of radius 
change. We then modify our earlier analyzes to con- 
sider other sources of the variations. Our analysis 
is preceded by an explanation of the meaning of the 
"seismic" radius inferred from f-modes. We then com- 
pare our results to those of Antia et al., and interpret 
the inferred rates in terms of magnetic field and tem- 
perature changes. 

Our interpretation of p-mode frequency changes is 
predicated on the work of Goldreich et al. (1991) who 
analyzed BBSO measurements from the rise of the 
previous cycle (cycle 22). We use our inference on the 



behavior of the sub-photospheric layers to constrain 
radius changes arising there. 

2. Frequency data from SOHO/MDI 

In the present study, we use 19 MDI data sets 
containing centroid frequencies determined from mea- 
surements made between May 1, 1996 and June 21, 
2000, with a break between June 16 and October 22, 
1998, when there was no contact with SOHO. The 
sets are typically 72-days long, except those immedi- 
ately before and after the break, which are shorter. 
The centroid frequencies, V£ in , were determined by 
the method described by Schou (1999). 

The sets contain between 112 and 203 f-mode fre- 
quencies, with earlier sets having more data. The 
maximum lvalue is 300 and the minima range from 
89 to 137. The number of p-mode frequencies range 
from 1589 to 1906. Again, the earlier sets are mode 
abundant. The p-modes range between £=0 and 200. 
The differences in mode composition are not impor- 
tant in the case of p-modes. The number of overlap- 
ping modes is large enough for a detailed study of 
frequency changes. In the case of f-modes, the differ- 
ence in the f-range may be important, and therefore 
in our study of the solar radius changes, we used only 
modes with I > 137. 

3. Inferences from f-mode frequency changes 
3.1. Helioseismic radius 

All helioseismic determinations of the solar radius 
to date have relied on the following asymptotic rela- 
tion for f-modes frequencies, 

u e 2 R ' U 

Antia et al. (2000) pointed out that using this rela- 
tion for modes with £ extending up to 300, as Dziem- 
bowski et al. (1998) did, is not justified because of 
significant departures from v oc R~ 15 are present in 
higher ^'s. The departure increases with £, which 
as Brown (1984) first suggested could be accounted 
for as an effect of turbulence in the upper convec- 
tive zone. Detailed models of this effect have been 
developed by Murawski & Roberts (1993a, 1993b). 
(For the most recent work on the subject see Medrck 
& Murawski 2000). However, surface magnetic fields 
may also have significant effects on f-mode frequencies 
(Evans & Roberts, 1990; Jain & Roberts). With these 
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two sources of perturbation to f-mode frequencies, we 
must contemplate solar cycle changes beyond that of 
a simple radius change. The relative contribution of 
the near-surface changes are expected to increase with 
I, because such changes should be inversely propor- 
tional to mode inertia, It, which sharply decreases 
with I. 

There is another problem in applying Eq.(l) in a 
search for the radius variations correlated with activ- 
ity. This problem follows from the fact that the in- 
duced modifications are quite non-uniform, and each 
f-mode has it is own radius, Re, which is given by 



Rl 



-1/3 



(2) 



With this definition, we get from the variational prin- 
ciple for oscillation frequencies (see Appendix) 



1 ,r ^ GM 



(3) 



where L — \J 1(1 + 1). This is a very accurate ex- 
pression. The relative departures from equality range 
from 2 x 1(T 4 at I = 100 to 8 x 10~ 5 at I = 300. In a 
linear approximation in terms of (R — Re), this is the 
same as the formula obtained by Gough (1993). 

For high degree modes, the f-mode radii are close 
to the solar radius. The values of Re/R range from 
0.9883 at I = 100 to 0.9946 at t = 300. While we 
have Re~R,& corresponding approximation for ARe 
is quite problematic. When the f-mode frequencies 
were used to refine the value of the radius for model- 
ing the sun, we could expect an approximate, homolo- 
gous relation, Re oc R. But such a relation cannot be 
expected in the case of the activity induced changes, 
which we believe to be confined to the outermost part 
of the sun. If the data show that Aug <xve, as Antia 
et al.(2000) found, then the simplest interpretation 
is that, indeed we have ARe being constant over the 
range of considered I- values. Then, the inferred value 
of AR in Eq.(l) refers to the range of depths beneath 
the photosphere corresponding to the range of £'s in 
the data sets. Antia et al. (2000) used modes in the 
100-200 range, which translates to 10-6 Mm in depth. 
Their finding implies that this layer was moved down- 
ward by about 5 km during the two years they con- 
sidered. The truth is, with these data, we cannot say 
anything about what happened in the layers above. 
Thus, we have no information about the evolution of 
the photosphcric radius of the sun. 



3.2. Formal determination of the rate of shrink- 
ing from f- modes 

To account for the effect of the near-surface changes 
on f-mode frequencies and possible differential changes, 
we modified Eq.(l) into 



ZARf Aj f 



(4) 



where ARf denotes the radius change inferred from 
a particular set of f- modes. For the calculation of It, 
we adopted the following normalization of the eigen- 
functions 

(£ r ) ph = 2 -^Ye m , (5) 



where £ r denotes radial displacement of the fluid ele- 
ment. With such a normalization, the values of It(v) 
are of the order of unity for p-modes in the 2-4 mtlz 
range. For the f-modes, the Ie values are significantly 
larger (see Table 1). 

We determined ARf and 7/ by a least-squares fit- 
ting. In Fig. 1, we show values of Ave for selected 
data sets. The first two sets were obtained at solar 
minimum. The 1999.4 set was taken near the mid- 
dle of the phase of rising activity, and the last is at 
the current maximum. Here, A denotes the difference 
between the solar data and the reference model. The 
reference solar model is that of Christensen-Dalsgaard 
et al.(1996). It has the same radius used by Dziem- 
bowski et al. (1998). That is, the model uses the 
"old", too large value of the solar radius (not that 
of Brown and Christensen-Dalsgaard 1998), and this 
is why the frequency differences are rather large. A 
small difference in the reference model radius is in- 
consequential for the inferred temporal changes. 

We assumed that A7/ is ^-independent, and as we 
see in Fig. 1, the solid line is a good fit to the data. 
The x 2 's vary from 1.2 (2000.4 set) to 1.84 (1996.4), 
except for the significantly worse fit (x 2 — 3.55) found 
for the 1998.9 set, which was the first one taken after 
the recovery of SOHO. We also tried fitting 7/ as 
a low-order polynomial depending on frequency, but 
this did not improve the fit. 

We see the departure from the linear relationship 
implied by the the radius adjustment sharply in- 
creases with v. Antia et al. (2000) considered only 
modes with v < 1.44 mHz, and it seems that the de- 
parture from a straight line is still small. However, 
this is somewhat misleading because we used a model 
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with much too large a radius. As we shall see in Ta- 
ble 1, at the level of changes of a few nanohertz (i.e. 
radius changes of a few km), the difference is quite 
significant. We emphasize that high i- modes arc im- 
portant because with increasing £, Rg approaches the 
solar radius. For such modes, including 7/ is essential, 
which implies that we have to rely on Eq.(4) rather 
than Eq.(l). With Eq.(l), we get a much poorer fit 
(x 2 = 4.4— 16.5) and the correction to the solar radius 
is larger by some 20 km. This illustrates the trade-off 
- increasing I moves us closer to the surface, but such 
high £'s are more strongly contaminated. 

In Fig. 2, we show the variations of the solar ra- 
dius and 7/ inferred from f-modes from the truncated 
data sets. The rise of the current activity cycle be- 
gan in 1997.4 which was marked by a sharp rise of the 
seismic activity indicators (Dzicmbowski ct al. 1998). 
A corresponding sharp rise of p-mode frequencies be- 
ginning at this time may be seen in Fig. 3 here. That 
is why we choose 1997.4 to begin our linear fits. We 
have no explanation, as yet, for the relatively large 
fluctuations in ARf which appear to have a one- year 
period. For comparison, we also show the result ob- 
tained when the 7^-term is ignored. There is a dif- 
ference, but not as large as one might anticipate by 
looking at Fig.l. The rate of radius decrease is only 
insignificantly higher than in our standard version, 
and the error is larger. 

In detail, we found from our linear fit, with the 7/, 



dR 
~dt 



(-1.51 ±0.31) km/y, 
and without the 7/-term, 



(6) 



dR _ 
~dt~ ( -' 



■1.82 ±0.64) km/y. 



The values are similar to those found by Antia et al. 
(2000). To make a closer comparison, we truncated 



Table 1: Contributions to f-mode frequency shifts 
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Fig. 1. — Differences between measured and calcu- 
lated f-mode frequencies. The error bars show esti- 
mated standard deviations of measured values. The 
dates correspond to the center of the individual 72- 
day long measurement periods. The solar model was 
calculated assuming i? Q = 695.991 Mm. The solid 
line represents the fit to Eq. (4) . The dashed straight 
line represents the part attributed to the difference 
between the solar radius and that adopted in the 
model. 
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DATE 

Fig. 2. — Upper panel: Variation of solar radius be- 
tween 1996.4 and 2000.4 inferred from f-mode fre- 
quencies with and without the 7/-term. Two straight 
lines represent linear fits to the data starting from 
1997.4 when the rise of cycle 23 began. Lower 
panel: Corresponding variation of 7/, which describes 
remaining near-surface contribution to f-mode fre- 
quency variations. 
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our data sets at £ = 200, and then we found 
dR 

— = (-1.80 ±0.38) km/y. 

Having in mind that we still miss modes between 
I = 100 and 137, it is fair to say that there is no dis- 
agreement between our findings and theirs, implying 
that at a depth of from 6 to 10 Mm the sun shrank by 
some 4 to 6 km during the rising phase of this activity 
cycle. 

How reliable is this finding? The main concern 
is the role of the near-surface perturbation and the 
cross-talk between the two terms on the right hand 
side of Eq.(4). In the lower panel of Fig. 2, we show 
the 7's. The linear fit for 7, which is visibly poorer, 
yields 

<% 
dt 

The relative contribution of the two terms to overall 
f-modc frequency variations depends on L In Table 1, 
we compare these two contributions, denoted byAfR 
and Av 7 for selected I- values. The increasing role of 
Av 1 is a consequence of decreasing mode inertia. It 
should be noted that Av 7 yields an appreciable contri- 
bution to Av even for modes with I < 200. Caution is 
necessary, but we will proceed further assuming that 
the effect is indeed real. 

3.3. Accounting for the rate of shrinkage 

Even as small as it seems, a shrinking of the sun's 
radius during the rising phase of activity is not easy 
to explain. To investigate, we write the Lagrangian 
change of the local radius in the form 



where 



(0.180 ±0.051) pHz/y. 



(7) 



Ar(r ) 



-I 



Ap ( x 
, P \ro 



dx, (8) 



where r& is the radius at the bottom of the layer per- 
turbed by activity, and r$ is the radius at a specified 
fractional mass, M r /M, at activity minimum and Ap 
denotes the horizontally averaged change of density. 
We obtain a more revealing form of Eq.(8) by express- 
ing Ap in terms of the averaged entropy and magnetic 
field changes. 

For the horizontally averaged gas pressure in the 
presence of a random magnetic field we have, after 
Goldreich et al. (1991), 



AP n 



-A(PPm), 



(9) 



Bl + Bl 



8tt 



is magnetic pressure and 

0=-^ 



Bl 



87rP m 

is a measure of the statistical anisotropy of the field. 

With the use of thermodynamical relations, we de- 
termine 

2 



Ar 



/ 

Jr h 



1 A(pP m ) AS 

5 + (-Pr) — 

a c p 



Fi 



dx, 



(10) 

where px denotes the logarithmic derivative of density 
at constant pressure. The remaining thermodynami- 
cal quantities have their standard meanings. At the 
relevant depths, the gas is nearly ideal. Thus, we may 
use pr = —1, 1/Ti = 0.6, and find 

AS 

Cr, 



f^-0.4^. 



The irradiance from an active sun is higher than av- 
erage. If the same is true about luminosity then we 
should have AS* < 0. Hence, a negative contribution 
to Ar. However, this must be very small. If AS refers 
to the whole convective zone then a 10 -3 luminosity 
increase translates to an annual decrease in AS/cp 
of 10~ 7 . Another possibility is an increase in the su- 
peradiabatic gradient, V con — V a d, but this seems un- 
likely too. The annual decrease of ARf = 1.5 km 
refers to the layer of r/R = 0.988 - 0.995. Thus, 
ARf must arise mostly beneath r = R137 = 0.988i?. 
At this depth, according to a mixing-length model, 
V con — V a d ~2x 10~ 4 , which rapidly decreases go- 
ing inward. We would need an order of magnitude 
increase in the superadiabatic gradient to account for 
our rate of shrinking. 

A more acceptable explanation would be a vari- 
ation in the magnetic field. The consequences of a 
magnetic field increase depend on 0. For a purely 
radial field (0 = —1), the increase implies contrac- 
tion. For an isotropic field (0 — 1/3) the increase 
implies expansion. The field geometry implying the 
minimum increase to account for the rate of the 
shrinking corresponds to = —1. Then, we have 



A < B > rms = y A(B^), and assuming a constant 
rate across the lower convective zone, we infer 



d< B> r 
df~ 



7.2kG/y. 
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The value at R137 may be reduced, for instance, to 1 
kG/y if one allows an exponential increase of the rate 
to about 43 kG/y at the base of the convection zone. 

Thus, what we have inferred from the f-mode fre- 
quency change is a non-trivial constraint on the in- 
ternal magnetic field change. Let us note that if 
the field increase were predominantly isotropic, we 
would see an expansion rather than a contraction. 
The field increase inferred from the residual (after 
removing the near surface contribution) part of the 
p-mode frequency change was about 60 kG at 25-100 
Mm (Dzicmbowski et al., 2000). This high value could 
be consistent with the shrinking rate only if /3 is close 
to zero, that is if the field is essentially force- free, 
which is not a likely possibility. Thus, we are now 
skeptical about the reality of that large field change 
we reported earlier. 

Our inference regarding the solar radius change 
is limited by the lack of accurate information about 
what happened in the outer 4 Mm of the solar interior. 
This is the region where we may expect the largest ac- 
tivity induced variations for two reasons. First, the 
rapid decline of gas pressure and second, the thermal 
structure of this layer is more susceptible to changes 
in the efficiency of the convective energy transport in- 
duced by the field changes. The f-mode data we have 
at hand provide some information about changes in 
this layer through the 7/. Similar, but much more 
accurate information is available in the p-mode data, 
which we now consider. 

4. Inference from p-mode frequency changes 




1996.6 1997.4 1998.2 1999.0 1999.8 



DATE 

Fig. 3. — Variation of the mean value of 7 with two 
versions of its polynomial dependence inferred from p- 
mode frequencies. As in the case of f-modes shown in 
Fig. 2, the linear fit corresponds to the data starting 
from 1997.4 when the rise of cycle 23 began. The 
error bars would be within the symbols. 



4.1. The near-surface source of the p-mode 
frequency changes 

The p-mode spectrum of MDI frequency data is 
about 13 times richer than that for f-modes. Unfor- 
tunately, p-modes are not directly useful for deter- 
mining changes in the solar radius. The simple rela- 
tion, v oc R~ 15 would be valid for p-modes only if 
the changes were homologous throughout the whole 
sun. This is far from the truth for the changes we 
are considering here. However, from p-modes one 
may make a much more precise determination, than 
from f-modes, of the near-surface perturbation. For 
p-modes, we call it Aj p , and it describes frequency 
changes caused by a variable perturbation localized 
near the surface. In the present application, however, 
taking into account the v dependence is required for 
an accurate fit. We express the dependence in the 
form of a Legendre polynomial series with argument 

_ v - [yi + v h )/2 

s — > 
Vh - vi 

where v\ and Vh denote the lowest and the highest 
frequencies in the data set. Thus, we write 

1 3 

Here A is with respect to the 1996.4 data set. The 
number J was increased until 7 Pj0 stabilized within 
the errors and x 2 stabilized. This occurred for J = 2. 
In Fig. 3, we plot A7 Pi o for J = and 2. Variations of 
7 P are indeed much more accurately determined than 
those of 7/. For J = 2, we find the rate 

= (0.149 ±0.008) yuHz/y. (12) 

The dependence of 7(f) yields an important con- 
straint on the localization of the source of solar cycle 
variations in p-mode frequencies. 

Following Goldrcich et al. (1991), we link the fre- 
quency change to the change of the mean squared 
magnetic field and a Lagrangian change of a single 
thermodynamic parameter. For the latter, we pre- 
fer to use temperature rather than entropy which was 
used by Goldreich et al. From Eqs. (14) and (15) of 
Goldreich et al., we get the following expression for 
the change of 7 P , 

1 f -> AT 
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+[i + ri(r P + i> P ) 
-/3(ri-i + r lPP )]AP ro }. 



(13) 



Here, we denote by Tp and T p , the logarithmic deriva- 
tives of Ti. The ideal gas equation cannot be used 
in the layers where most of the contribution to Ajp 
arises. 

Goldreich ct al. (1991) explained the p-mode fre- 
quency changes during the rising phase of cycle 22 
in terms of magnetic field and temperature changes, 
with former being dominant and causing the fre- 
quency increase. They invoked a chromospheric tem- 
perature increase to explain the reversal in the in- 
creasing trend in A(^). We do not see such a trend in 
our data. Thus, as a first guess we interpret A"f p in 
terms of magnetic field changes. Later, we will discuss 
other sources of the p-mode frequency changes. 

We considered two values of /?, -1 and 1/3, and 
the following form for the depth, D, dependence of 
magnetic field increase 



A <B> r 



B b 

B b + \( Db _ 
A <B> rms 



D-D b 



-D 

(A) 



if D > D b 
if D t < D < D b 
if D < D t 



where D m = — 0.485Mm denotes D at the tempera- 
ture minimum, and B b , D b , and A were determined 
by fitting the three terms in the series given by Eq. 
(11). For Dt we adopted either D m or 0. 

In Fig. 4 we show two examples of the field's chang- 
ing behavior that would be consistent with the ob- 
served 7's, and compare them with two cases that are 
clearly inconsistent. One of the two inconsistent cases 
is a depth independent increase, and the other is an 
example of the field gradually increasing to about 3 
kG at 8Mm. In all four examples, we used (3 = — 1. 
We see that indeed the 7 P (^) provides a strong con- 
straint on the localization of the source of frequency 
changes, but clearly not a unique answer. For the two 
fitted cases, the inferred values of B b are 290 and 250 
G. Corresponding values of B p ^ = A < B > rms (0) 
are 62 and 94 G. An equally good fit was obtained 
with the choice (3 = 1/3. Data on the three models 
of the magnetic field change fitting A-y p (v) data are 
given in Table 2. The result for (5 = 1/3 is not sig- 
nificantly different from that found by Goldreich et 
al. (1991). To explain the p-mode frequency increase 
between minimum and maximum, we require an in- 
crease of the rms magnetic field growing from 0.2 kG 
in the photosphere to 0.84 kG at 4.25 Mm. The cor- 
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Fig. 4. — In the top panel, points with the error bars 
represent A70 {y) inferred from p-mode frequency dif- 
ference between 2000.4 and 1997.4. The lines corre- 
spond to various distributions of the averaged mag- 
netic field, shown in the middle panel. The solid and 
dotted lines are within the error bars in the top panel. 
The bottom panel shows the relative temperature de- 
creases required to cause similar frequency shifts as 
the corresponding magnetic field increases 
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responding numbers of Goldreich et al. are 0.25 and 
1 kG. 

In the bottom panel of Fig. 4, we plot the rela- 
tive temperature changes, which gives the same local 
contributions to A^ p as the corresponding changes in 
the magnetic field. We see that the required change of 
temperature is unacceptably large in the atmospheric 
layers. However, in sub-photospheric layers we can- 
not exclude the rms AT/T at the 10 -3 level. Such a 
temperature decrease would be a significant contribu- 
tor to the observed frequency increase. Bruggen and 
Spruit (2000) argue that one expects a lower subsur- 
face temperature from an increasing magnetic field, 
and that the effect should be searched for by means 
of hclioscismology. A contribution from temperature 
decrease would lower the requirement for the mag- 
netic field increase in the sub-photospheric layers. 

Yet another potential contributor to the frequency 
increase is a decrease in the turbulent velocity. Roughly, 
the relative change in the turbulent velocity, Avt jvt — 
q, has the same effect as a relative temperature change 
AT/T = O.HqM 2 , where M is the turbulent Mach 
number. In the sub-photospheric layers, M. is in the 
0.1-1 range. Thus, the effect may be significant, and 
we may expect a decrease in v t , with increasing ac- 
tivity, because the magnetic field should inhibit con- 
vection. 

4.2. Shrinking or expanding of the outermost 
layers 

In Table 2, we provide the values of the contri- 
bution to the rate of the photosphcric radius change 
due to the magnetic field increase inferred from the 7 P 
changes. We emphasize that the rate does not refer 
to photosphere but to the mass point corresponding 
to the unperturbed (solar minimum) photosphere and 
that the value does not include the part that was in- 
ferred from f-mode frequency changes. 

The solar photosphere is defined as a surface of 
specified optical depth r p h = M p hR, where M p h is 



column-mass depth and R is the mean opacity in the 
atmosphere, or, which is closely related, the place 
where the local temperature equals the effective tem- 
perature. Thus, if we want to assess the rate of move- 
ment of the photosphere, we have to take into account 
a possible change in k. To keep r p h unchanged, an 
additional radius shift of — An(dr / dR) p h is needed. 
Hence, the rate of the photosphere's change may be 
assessed as 

dt ~{dt) ph U<W Ph ' [ ' 

An estimate shows that the second term may not be 
negligible if (AT/T) ph - 10" 3 . 

The connection between i? ph and the solar disk 
radius, Rd, determined from the inflection point in 
the limb-darkening function was discussed recently by 
Brown and Christensen-Dalsgaard (1998). They find 
Rd — Rph ~ 500 km. Again a 10~ 3 temperature per- 
turbation within the atmosphere may be significant at 
the level of the radius changes discussed here. Thus, 
the difference between the solar radius variations in- 
ferred by means of seismology and photometry has to 
be kept in mind when a detailed comparison is made. 

The total value of (dR/dt) p h may be estimated as 
the sum of -1.5 km/s inferred from the f-mode data 
and one of the values inferred from p-modes data 
shown in Table 2. These are model dependent. We 
do not expect that by including the nonmagnetic con- 
tributors to Aj p , we would infer rates significantly 
beyond the range of values quoted in this table. We 
note that the net effect may imply both contraction 
and expansion. Possible net values of (dR/dt)^ range 
from -3 to 1 km/y. 

Finally, we point out that there is no contradiction 
between our inferences from f- and p-mode frequency 
changes. The effect of the field increases needed to 
account for the dRf/dt value have a negligible effect 
on p-modes frequencies, if the outward decrease of 
d < B > rms /dt from the bottom of the convective 
zone is steep enough. 



Table 2: Inference from p-mode frequency changes 5. Conclusions 

between 1997.4 and 2000.4 . 

Results of our analysis of t-modc frequency confirm 

.the evidence, first found by Antia and Basu (2000), 

for a contraction of the sun's outer layers during the 

rising phase of the magnetic activity. The rate we 

. determine is 1.5 km/y and is only somewhat different 

than found by our predecessors. We pointed out, how- 





A 


D b 


D t 


B b [G] 


B p h[G] 


(dR/dt) ph 


-1 


0.575 


3.00 


-0.485 


29 


62 


-1.3 km/y 


-1 


1.15 


1.27 





25 


94 


-1.1 km/y 


1/3 


0.623 


4.25 


-0.485 


84 


200 


2.3 km/y 
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ever, that there may be another interpretation for the 
observed frequency variation. Further, we stressed 
that the rate does not refer to the surface radius, but 
to the layer at 4 - 8 Mm depth below the photosphere. 
In spite of the fact that the dispersion relation for 
high-degree f-modes approaches that for the surface 
gravity waves, the two types of modes are essentially 
different. While the latter are discontinuity modes 
which see the same gravity for each horizontal wave 
number, the f-modes see different effective gravities 
depending on I. 

The rate of shrinking is most easily explained as 
resulting from the rise of the radial component of 
the random magnetic field beneath a depth of 8 Mm. 
There is an integral constraint on the magnetic field 
which may be the most important finding from the 
data on f-mode frequency changes. To account for 
the shrinking rate, we need an increase in the radial 
component of the random magnetic field with a mod- 
est annual rate. An isotropic increase would imply 
an expansion in the f-mode region. We pointed out 
that this new constraint is likely to be in conflict with 
the much larger change of the interior field inferred 
by Dzicmbowski et al.( 2000) from the inversion of 
p-mode frequency changes. 

The p-mode frequency change may be accounted 
for in terms of magnetic field changes. Our analysis 
was based on the formalism of Goldreich et al. (1991), 
and we found similar implications regarding the re- 
quired field increase as these authors, who analyzed 
BBSO data from the previous solar maximum. In par- 
ticular, the increase must be larger below the photo- 
sphere than in the atmosphere, if this is the sole effect 
causing p-mode frequency changes. We pointed out, 
however, that a temperature decrease and/or decrease 
of turbulent velocity in sub-photosphcric layers could 
be significant contributors to the frequency decrease. 
Depending on the field anisotropy, the changes in the 
outermost layers may lead to additional shrinking or 
to net expansion. 

Our estimated rates of radius change during the 
rise of cycle 23 range from -3 to 1 km/y. This dif- 
fers from the rate of about 5.9±0.7 km/y determined 
by Emilio et al. (2000) from the direct radius mea- 
surements based on SOHO/MDI intensity data. Per- 
haps the difference may be explained by the differ- 
ence between dRd/dt and our (dR/ dt)^. Both re- 
sults, however, imply a negligible contribution of the 
radius change to the solar irradiance variations. Fur- 
thermore, the two estimates of the radius change be- 



tween maximum and minimum activity are by two 
orders of magnitude less than found by Noel (1997) 
from his measurements with the astrolabe of Santi- 
ago. He finds the difference between the 1991 (pre- 
vious maximum) and 1996 radii which is exceeding 
700 km. The data from the Solar Diameter Moni- 
tor (Brown & Christensen-Dalsgaard, 1998) are in- 
consistent with such large variations, although there 
is a hint of possible radius increase during 1987 of 
some 30-40 km. On the other hand, a theoretical con- 
straint on radius given by Spruit (1994) is even tighter 
than than that from hclioscismology. The number 
he quotes for the maximum to minimum difference is 
2 x 10- 7 i? Q = 0.14 km. 
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APPENDIX 

We assume the Cowling approximation and write the 
equation for adiabatic oscillation in the following form 

pu 2 $, = VP' + p'ge r = F£. (Al) 

The notation here is a standard one and docs not 
require explanation. The variational expression for 
eigenfrequencies is 

2 _ [dh^Fj = K 

10 jd^ P \^ ~ r {A2) 

The f-modes are nearly incompressible. Thus, for the 
approximate £ to be used in this expression, we as- 
sume 

V-£ = 0. (A3). 

Then for the Eulcrian perturbation of density and 
pressure, we have 

P' = gpt r and p' = (A4) 

We express in a standard way the displacement eigen- 
vector in terms of the spherical harmonics, 

£= [y(r)e r + z(r)V}Y e m . (A5) 
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With this expression (A3) becomes 



dr r 



L - = 



and the integrals in (A2) become 

fR 

2„2\ 



(y 2 + L 2 z 2 )pr 2 dr 



(A5) 



(A6) 



and 
K 



R r 



2L 2 yz 



ding 
dlnr 



y 



-pr 2 dr, (A7) 



where we made use of (A4). We may use 
dlnq 

± = _2 

dlnr 

because for modes considered here the logarithmic 
derivative of the local mass, M r , is less than 10 -2 
in the layers contributing to / and K, which implies 
less than a 10 -4 fractional contribution to frequencies. 
From the ratio of radial to horizontal component of 
(Al), we obtain approximately 



L 



V 



and, taking into account the inner boundary condi- 
tion, y = Lz. Now, we have from (A6) 



1=2 y 2 pr 2 dr 
Jo 



and from (A7) 



K = 2(L - 2) ( y 2 -pr 2 dr 
Jo r 



(A8) 



(A9). 



Eqs.(2) and (3) follow immediately from (A2), (A8), 
and (A9) upon setting M r = M and Li = 2-kv. 

REFERENCES 

Antia, H.M. 1998, A&A, 330, 336 

Antia, H.M., Basu, S., Pintar, J. & Pohl, B. 2000, 
Solar Phys., 192, 459 

Brown, T.M. 1984, Science, 226, 687 

Brown, T.M. & Christensen-Dalsgaard, J. 1998, ApJ, 
500, L195 

Briiggen M. fc Spruit H.C., 2000, Solar Phys., in 



Christensen-Dalsgaard, J. et al. 1996, Science, 272, 
1286 

Dziembowski, W. A., Goode, P. R., DiMauro, M.P., 
Kosovichev, A.G. & Schou, J. 1998, ApJ, 509, 456 

Dziembowski, W. A., Goode, P. R., Kosovichev, A.G. 
& Schou, J., 2000, ApJ, 537, 1026 

Emilio, M., Kuhn, J.R., Bush, R.I., & Scherrer, P. 
2000, ApJ, 543, 1007 

Evans, D.J. & Roberts, B. 1990, ApJ, 356, 704 

Goldreich,P., Murray, N., Willette, G. & Kumar, P. 
1991, ApJ, 370, 752 

Gough, D.O. 1994, in Astrophysical Fluid Dynamics., 
ed. J. -P. Zahn & J. Zinn-Justin (Amsterdam: El- 
sevier) p. 399 

Jain & Roberts, B. 1990, A&A, 286, 243 

Medrek, M. & Murawski, K. 2000, ApJ, 529, 548 

Murawski, K. & Roberts, B. 1993a, A&A, 272, 595 

- 1993b, A&A, 272, 601 

Noel, F. 1997, A&A, 325, 825 

Schou, J. 1999,ApJ, 523, L181 

Schou, J., Kosovichev, A. G., Goode, P. R. & Dziem- 
bowski, W. A. 1997, ApJ, 489, L197 

Spruit H.C. 1994 in The solar engine and its influence 
on terrestrial atmosphere and climate, NATO ASI 
Series I., vol.25, p.107 



press, astro-ph/0004140 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v4.0. 



10 



